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Abstract 

The diagonal hydrodynamic reductions of a hierarchy of integrable 
hydrodynamic chains are expUcitly characterized. Their compatibihty 
with previously introduced reductions of differential type is analyzed 
and their associated class of hodograph solutions is discussed. 
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1 Introduction 



In a series of papers PQ-^ we have considered an infinite hierarchy of inte- 
grable systems [Hj which admits many interesting (1 + l)-dimensional reduc- 
tions hke the Burgers, KdV and NLS hierarchies as well as different classes of 
energy- dependent hierarchies |E]-jn| including the Camassa-Holm model also 
[TUj . It motivated the term universal hierarchy we proposed in p. 

This universal hierarchy can be defined in terms of a generating function 
G = G{X, x) depending on an spectral parameter A and an infinite set of 
variables x :={... , x_i,Xo, xi, . . .), which admits expansions 

G = IH h H , A ^ oo, (1) 

G = boix)+bi{x)X + b2ix)X^ + --- , A^O. (2) 
The hierarchy is provided by the system of flows 

9„G=(A„,G), neZ, (3) 

where 

{U,V) := Uid,V)-id,U)V, 

and 

^n: = A" + (7i(a;)A"-i + --- + (7n-i(a;)A + (7„(a;), n>0 (4) 



A_.: = M^ + M^ + ... + ^, n>0. (5) 

A" A"-i A ^ ^ 

In terms of the coefficients {gn{x)}n>i U{^n(^)}n>o the system Q becomes 
a hierarchy of hydrodynamic chains. 

An alternative and useful formulation of Q is obtained by introducing 
the generating function 

H -.^ i (6) 

Thus Eq.(jHl) is equivalent to 

d^H = d,(^A^H^, neZ, (7) 
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which means, in particular, that the coefficients {hn{x)}„>i supply an infinite 
set of conservation laws for (jS)). 

The hierarchy Q forms a compatible system. Indeed, as a consequence 
of one derives the consistency conditions 

dnAm - dmAn = {An, Am), 12,171 EZ, (8) 

It is important to notice that (jHJ implies that the pencil of differential 
forms 

..(A) ■.= J2iHAn)dxn, (9) 

nGZ 

is closed with respect to the variables x since it satisfies 

dm{HAn)—dn{HAm) = d^H An — dnH A^ + H(dmAn — dnAm) 
=dx{HAm)An — dx{HAn)Am + H {An, Am) = 0. 

The potential function Q = Q{\, x) corresponding to u 

dQ = uj, (10) 
leads to another useful formulation of our hierarchy. Indeed, according to (jOj) 

dnQ = And,Q, neZ. (11) 
Moreover, (fTTj) is completely determined by Q as 

"-^^ ^''^ 

In Section 2 of this paper the formulation (|lip is applied to prove that 
includes multidimensional models such as 

^yy ^y^xz UxyUz, (-^^) 
Uzx '^x^yy ^xy'^y, ('^^) 

^) = m . (16) 
-) = f^) • (17) 

Ux't \Ur'x 
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Section 3 deals with the main aim of the present paper: to characterize the 
hydrodynamic reductions of 0- These reductions are given by the solutions 
of of the form G = G{X, R), where R = {R^, . . . , R^) denotes a finite 
set of functions [Riemann invariants) satisfying a system of hydrodynamic 
equations of diagonal form 

dnR' = AUR)d^R\ neZ. (18) 

This type of reductions appeared in the context of the dispersionless KP 
hierarchy PI]-^^ and has been also used in to characterize the integra- 
bility of (2 + l)-dimensional quasilinear systems. In the present paper we 
study these reductions for the whole set of flows of the hierarchy (jH]) and we 
obtain their explicit form . Our results are summarized in Theorem 1. It 
should be noticed that a similar analysis for the first member (ti-fiow) of 

has been recently performed in ^B]. The compatibility between hydro- 
dynamic and differential reductions considered in part 2 of Section 3. The 
paper finishes with a discussion of the solutions supplied by the generalized 
hodograph method 17]. 

The following notation conventions are henceforth used. Firstly, Goo and 
Go stand for the expansions (0) Q of G as A ^ oo and A ^ , respectively. 
Furthermore, let V be the space of formal Laurent series 

oo 

= ^ a„A". 

n=— oo 

We will denote by Vr,s {r < s) the subspaces of elements 

s 

F = 5^a„A^ 

n=r 

and by P^.s : V ^ Yr,s the corresponding projectors. Given 1/ G V we will 
also denote 

In particular, notice that we can write 

An = (a"Goo) , A_n = (a-"Go) , n > 1. 

V / 0,+oo V / -oo,-l 
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2 Integrable models arising in the hierarchy 
2.1 Multidimensional models 

If we use the potential function Q{X, x) for the differential form (jH)) 

dQ = Y,iHA,,)dxn, (19) 

then from one readily deduces that Q admits expansions of the form 

<y = 2^Ax„H h H , A ^ oo, (20) 

n>0 



^ = E5+Po(^)+Pi(^)^ + P2NA2 + ---, A^O. (21) 

n>l 

By substituting into (fT^ and by identifying coefficients of equal 
powers of A one obtains formulas for the differentials of the functions qn{x) 
and Pn{x) in terms of the coefficients of the expansions (jT)) and 0- For 
example, the simplest ones are 

dqi= 6odx_i + ^(&„ dx_„_i - dx^-i), (22) 

n>l 

dpo = -^(dxo + ^ign dxn-bn dx^n))- (23) 

n>l 

They imply 

dqi = ^^(^da;_i - d^+iPo dxn^ (24) 



dpo = f dxp - y^(9n-igi dxnV (25) 

Permutability of crossing derivatives of qi and po in these identities lead at 
once to multidimensional nonlinear equations for the functions Qn and 6„. 
For example, starting from 

dmd-iqi = d^idmqi, m ^ -1, 
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and using (j24j) . the following nonlinear equation results 

dndoPo = d^po {d^idn+m) - {d-Apo) D^+iPo, n ^ -1. (26) 
In the same way, from the crossing relation 

dmdnqi = dndmqi, 171,12 -1, 

and ()24|l we get the following nonlinear equation 

3j^)=aJ^), m,n^-l. (27) 

The same type of equations can be derived for qi. The different choices 
available for n, m in the equations and ^Fi\ give rise to the models 
(113- (HZI). 

2.2 2-dimensional integrable models 

In PP-|S] we developed a theory of differential reductions of our hierarchy 
based on imposing differential constraints on G ^ Goo of the form 

(^(A,G,G,,C,...)) =0, x:=xo. (28) 

In particular the following three classes of reductions associated to arbitrary 
polynomials a = a(A) in A were characterized: 

Zero-order reductions 

a{X)G = U{X,x), U:=(a{X)G] , (29) 

V / 0, + OD 

First-order reductions 

G, + a{X)=U{X,x)G, U:=(^) , (30) 

VG/ 0,+oo 
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Second-order reductions 



^GG^,-^Gl + a{\) = Ui\,x)G', U:=[-^)^^^. (31) 

The first-order reduction for a linear function a(A) determines the Burgers 
hierarchy. On the other hand, under the differential constraints (jHlll the 
hierarchy ()18|) describes the KdV hierarchy and its generalizations associated 
to energy- dependent Schrodinger spectral problems. In particular, the linear 
and quadratic choices for a(A) lead to the KdV (Korteweg-deVries) and NLS 
(Nonlinear-Schrodinger) hierarchies, respectively. Indeed, if we define the 
functions ip{X, x) by 

ij{X,x) := exp(D;V), := ± H, (32) 

then from ()3H) it is straightforward to deduce that 

= An{~D^logH ± ^/aH)^jJ - ^A„,^^/' 

^xx = (0. + 0')^ = {{D,, H} + aH')^ = U^. 

In other words, the functions ip are wave functions for the integrable hierar- 
chies associated to energy-dependent Schrodinger problems. The evolution 
law of the potential function U under the flows |TH|) can be determined from 
the equation 

= —A-n^xxx "I" 2?7 An^x ~^ UxA^, (^3) 

which arises as an straightforward consequence of (fT5|) and (jHH). 

Additional reduced hierarchies including nonlinear integrable models such 
as the Camassa-Holm equation can be also deduced (see |lj). 
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3 Hydrodynamic reductions and solutions 
3.1 Hydrodynamic reductions 

Let us consider now the hydrodynamic reductions of Q. We look for classes 
a solutions 

G = G{X,R), (34) 

of where R = {R^, . . . , R^) satisfies a infinite system of hydrodynamic 
equations of diagonal form (jl8j) . Our aim is to characterize both the form of 
G = G{X, R) and the characteristic speeds defining the system (fTHj) . By 
substituting (jH^j) into Q then by using (fTHj) the identification of coefficients 
of the derivatives dxR^, i = 1, . . . , N implies 

(AG)A^ = v4„(AG') - (AA^)^, 1 < z < A^, neZ (35) 

where 



In addition to these equations we impose the requirement of the commuta- 
tivity of the flows (jl8|) . which is equivalent to the following restrictions on 
the characteristic speeds 



DA* DA* 



M - A*„ Ha - A* 



2 7^ j, m ^ n. (36) 



We start our analysis by considering the positive flows n>l. From ()35|] 
we get the following system for G ~ Goo 

(AG)At, = (A"G)_oo,o(AG) - (A(A"G)_oo,o)G', n>l. (37) 

By substituting in these equations the expansion G ~ Goo and identifying 
the coefficients in ^ and we get 

A^7n+i = A^A^7i, (38) 
Dign+2 = A^A5'2 + fi'„+i(A5'i) - iD^gn+l)9l■ (39) 

As an immediate consequence it follows 

A^+i = (7„+i+A^(Ai-(7i), (40) 
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which imphes 

Ai = AjX = A\-gA, n>l. (41) 



We now look for a system characterizing A^^^ and gi. To this end we 
differentiate PHj) and find 

= DiDjgiAi + D.g^DiAi^, 



so that 



D.,g, = -^^D^g, + -^^D,g,, 
Ai -Al Al- Ai 



n 



and from (jHUj) we may write 



A{-A\ Al-A{ 



D,,g,= !, ' \, Dj9i, J- (42) 



On the other hand according to (j4(jp 

Ai, = g2 + A{{A\-gi), 
so that from and with the help of ()H8j) we find 

DjA\ _ D,Ai _D,A\{2A\-g,) + {A{-A\)D,g, 



« 7^ J, 



A{ - A\ Ai - A| (Ai - Al) (Ai + A^ - g,) 

which reduces to 

DjA\ = D^g^, i ^ J. 

In this way 

A\ = g, + f{R'), (43) 
where the functions are arbitrary. By using this result in ()42j) it follows 

Dijgi =0, i^j, 

so that 

TV 

g, = J2h'{R'), (44) 

k=l 

where the functions h'^ are arbitrary. 
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By using these results we may determine G{X, R) since from the equation 
dnni) with n = l 

where /i* := Dih^. Thus we get 

G{X, R) = exp ( f f J^^^ d B^) , (45) 

where the undefinite integrations are determined up to a function of A decay- 
ing at A — > oo. The expression fl45p coincides with the generating function 
of the conservation laws densities found in |16j . 

Let us now determine the characteristic speeds for the negative flows 
n < —1. The equations ()H5j] imply the following system for G ^ Go 

{DiG)K\ = (A-"G)o,oc(AG) - (A(A-"G)o,oc)G, n > 1. (46) 

Then by inserting the expansion G Gq of © and identifying the coefficients 
in A° and A we get 

A&n = (&n + ALjAln6o, (47) 

6oAfen+l = (fen + K)D^hl + 6n+l(A&o) " {D.h^)^. (48) 

It implies the following recurrence relation for the characteristic speeds 

AVi = ^(&n + AL„), (49) 

Oo 

which implies 

AL„ = A_„(a = ^), n>l. (50) 
The only unknown now is K^__i since according to ()45j) 

To find A^i we use ()49|1 for n = 1 

A-2 = ^(&i+A-i), 
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and the commutativity condition for the n = — 1 and n = —2 flows 

Aii-A*li ~ Ai^-ALs' 
Thus by ehminating AL2 one finds at once 

Dj In Ail = Dj In 60, j- 

Therefore, 

Al, = bog\R'), (51) 

with being arbitrary functions. 

At this point we have determined all the unknowns of our problem. How- 
ever, in our calculation we only used a subset of the equations required, so we 
must prove that our solution satisfies the full system of equations (j^ - lfHU)) . 

Let us begin with (jH^ for n > 1 

AlnG= n>l, 

/in 

which according to (HH), (jlHjl - ipSj) reads 

h 

A-/ 

To prove this identity we express v4„(A) as 



AAn(A) = (a„(A) - An{X = f )) . (52) 



A„(A) = (A"G)o,oo = / ^^3^dA = -i./ ^C^d~X 

2vr2L A- A 27r«L A-A 



i=l 



where 700 is a positively oriented closed loop around 00 in the complex plane 
of A (A and /* are assumed to lie inside the loop). Now by differentiating 
this expression with respect to i?* one finds 

D.MX)^±[ J^^5^,X^±[ .^^M^) dA 
27rz J^^ (A - A) (A - f ) Svrz J^^ (A - A) (A - f ) 



^ 'A„(A)-A„(A = f) 



x-f 
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which proves This identity leads also at once to since it implies 



An{X = p)-A4X = f) p-f' 
which means that 

- n> 1. 



Ai - Ai Ai - A\ 
Let us consider now (jH^j) for n < — 1 

Di\nG= ,^'^7. , n>l. 
From (Uni), dSni) and (jSIl) it takes the form 

D.A^niX) = (a_„(A) - A_„(A = ^)) . (53) 

In order to proof ()53|) we express A_n{X) in the form 

271-2 A - A 27r2 A - A 

where 70 is a closed small loop with negative orientation around A = (A 
and /* are assumed to lie outside the loop ). By differentiating with respect 
to i?* it yields 

AA-n(A) = — / ^^^dA 



27Ti 



70 



(A-A)(A-f) 



^^^-"^^^ dA = -A_(A.„(A)-A_„(A = r) 



2vr2 7^„ (A-A)(A-f) A-/ 
Hence (jSSl) holds if we set 
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The remaining commutativity conditions (j36p follow at once. 

Therefore we may summarize our analysis in the next theorem 

Theorem 1. The hydrodynamic reductions of the hierarchy (jS)) are deter- 
mined by 

G{X, R) = exp ( f f -^^^^^ d r) , (54) 



dr^R' = KiR) d.R\ A'r.iR) ■■= A„(A = f{R')), (55) 

where the functions h^ and /* are arbitrary. 

Since the systems (|18|) are invariant under local transformations of the 
form i?* — i> R^{R^), without loss of generality we will henceforth set 

h\R') = R\ (56) 

so that the form of the reduced generating function is 

A^TW)^' ^^^^ 

In this way some of the simplest hydrodynamic reductions ()54|) are given 



G(A, i?) = exp ( ^ / 
i=i ^ 



by 



rm = c,, G = exp(^^-— ), (58) 

i=l * 

G^fl{^\ (59) 

1=1 

We also notice that the characteristic velocities of the hydrodynamic sys- 
tems ()55|) can be written in terms of the Schur polynomials 



n>l n>0 
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as 

n 
j=0 

n-l . 

AL„(i?)=exp(-/o)5^5,(/_i,...,/_,)(r(i?0)"' n>0, 

j=0 

where 

1=1 

3.2 Compatibility with differential reductions 

A natural question is to find the hydro dynamic reductions compatible with 
the differential reductions ()30|) and ()31|) . Let us prove the following result: 

Theorem 2. The only hydrodynamic reductions ()57|) compatible with either 
first or second order differential constraints are characterized by 

r{R) = -R + c,, z = l,...,Ar, 

which correspond to generating functions of the form 

N 

G{\R)=a{\)\{{\ + R -a). 

i=l 

Proof. If we substitute (|Ffj) into the differential constraint for second- 
order differential reductions we get 

a(A)_ d^^R' \ 1 d.^R' \2 

~ ~2^V (A - f )2 ^ J ^^'^ ' 

This means that G has a simple zero at each A = /*(-R*) so that from (jKTj) 
we have 

where H is different from zero at A = f^{R^). If we now differentiate with 
respect to i?* we deduce that 
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Therefore the statement of the theorem for second-order differential con- 
straints follows. The corresponding proof for first-order constraints is simi- 
lar. □ 



3.3 Hodograph solutions 

The general solution of the infinite system is provided by the implicit 
generalized hodograph formula [T7j 

x+ Yl KiR) ^n = r{R), 2 = i,...,iv, (60) 

nez-{o} 

where the functions F* are the general solution of the linear system 

D r* 1 
' - ^ + 2- (61) 



_ pi ji _ ji 
By introducing the potential function $(-R) 

r = ^ = i,...,iv, 

the system reduces to the Laplace type form 

(f - f )A/^,$ = A<f - 1 + 3. (62) 

the general solution of which depends on arbitrary functions of one vari- 
able. 

In particular it is immediate to deduce ^U] that the generating function 
()57|) provides a one-parameter family of solutions of ()62|) . Thus we can 
produce important solutions of from linear superpositions of G{\,R). 
For example 

Aj, = A$, $:=— / A"G'(A,i2)dA, n>0, 



2m 



7oo 



AL„ = A$, $ := — / A^"G(A, R)dX, n>0. 
2'Ki 
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Furthermore, in several important cases the general solution of ()62|1 can be 
written in terms of G{X, R). 

Examples 
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If we set 

nR) = -, n = l,2,..., 

n 

then the generating function (non normahzed as A ^ oo) is 

N 

G(A,i2) = n(A 



n 

1=1 



Hence, by taking for each i = 1, . . . , a closed loop 7, in the complex A- 
plane with positive orientation around Aj = the general solution of 
can be expressed as 



N 

$:=— ^ / 0,(A)G(A,i?)dA, 

i=l >^7. 



where the functions 0j(A) are arbitrary. For example, if n = 1 it takes the 
form 

^ 1 

1=1 k^i 

Similarly we may deal with the case 

f{R') = --, n=l,2,..., 
n 

which leads to 



G(A,i?) = n(A + 



n 

i=l 



Now to generate the general solution of ()62|) we take for each i = 1, . . . ,N 
a path ■ji{R'^) in the complex A-plane ending at Aj = —R^/n so that we can 
write 



1 ^ r 

$:=— V / 0,(A)G(A,i?)dA, 



where the functions 0i(A) are arbitrary. 
Let us consider in detail the case 



r{R) = -R\ 2 = 1,2. 
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One finds 



A, 



with 6i{X) being arbitrary functions. A normalized generating function is 
given by 

G{X,R) = (^x + x^y ' ^otO. (63) 

We may characterize the general solution of the hydrodynamic flows (I55|) 
corresponding to n = — 1 and n = —2 by means of the hodograph formula 

x + A_i{X = -R')y + A_2{X = ~R')z = Di^, z = l,2, (64) 

where y := z := X-2. From (jH!^ one calculates 

6o = c^R^R^, hi = c\R^ + R^) - 2c^R'R\ c := ^, 

Ao 

and gets that the system reads 

(x - c'z) - c\y - 2cz)R^ = ei{~R^){R^ - R^) - ^i(-i?^) - ^2(-i?'), 

(65) 

(x - c^z) - c\y - 2cz)R' = ^2(-i?')(i?' - R^) - Oi{-R^) - O^i^-R^). 
In particular, it implies 

R = R{x — (?z, y — 2cz), 
so that R is constant on the straight lines 

X := (x, y, z) = (xo, yo, 0) + (c^ + 2c, l)s. 
We notice that starting from these solutions we may generate solutions 



u{x,y,z) 
of the nonlinear equation 



j\o{R)dy + j h{R)dz, 
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